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<^ Abstract 

C^ We report on the non-trivial nanoscale kinetics of the deposition of novel complex oxide heterostructures 

composed of a unit-cell thick correlated metal LaNiOs and dielectric LaAlOa- The multilayers demonstrate 

^^ exceptionally good crystallinity and surface morphology maintained over the large number of layers, as 

^ confirmed by AFM, RHEED, and synchrotron X-ray diffraction. To elucidate the physics behind the growth, 

H 

"^ the temperature of the substrate and the deposition rate were varied over a wide range and the results were 

Cd treated in the framework of a two-layer model. These results are of fundamental importance for synthesis 
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of new phases of complex oxide heterostructures. 



Recently, novel many-body phenomena emerging at the oxide-oxide interface have attracted 
considerable attention motivated by the intriguing possibility that such interfacial phenomena can 
form the basis for new materials with greatly enhanced functionality [[1]-S1. While the multiple and 
coupled degrees of freedom (lattice, spin, charge and orbital) potentially offer a means for manipu- 
lating and controlling novel functionalities, an understanding of the mechanisms of layer-by-layer 
(LBL) growth of such heterostructures (HSs) with exquisite tunability and quality control still 
present a formidable challenge. Naturally the challenge of finding parameters for LBL growth in- 
creases manyfold, when the multilayers are subjected to large strain[:5], broken translation symme- 
tries of the interface, and dimensionality confined to a single unit cell[[6l. During the past several 
years, extensive work on pulsed laser interval deposition combined with high pressure reflection 
high energy electron diffraction (HP RHEED) opened up new prospects for the stabilization of 
unusual material phases [|3. 

Growing unit-cell thin complex oxide HSs in planar geometry is particularly challenging task. 
In this letter, we report on LBL growth of novel HSs having unit-cell-thick (uc) repetitions of 
correlated metal LaNiOs (LNO) interlaid with layers of the dielectric LaAlOs (LAO). The obtained 
results demonstrate that the growth kinetics of such luc LNO/luc LAO superlattice (SL) has 
several interesting features leading to important fundamental conclusions applicable to growth 
of a wider range of complex-oxide HSs. Specifically, our findings elucidate that, in addition 
to temperature, modulated flux can be employed as a powerful kinetic handle for the fine-tuned 
manipulation of sub-monolayer growth regimes. 

The LNO/LAO SLs have been grown on atomically flat TiOz-terminated SrTiO3(001) (STO) 
single-crystal substrates prepared by the recently developed wet-etch procedure[[8l to minimize 
surface and near-surface electronic defects. The PLD system was equipped with a newly developed 
high pressure RHEED system operating in the background O2 pressure of up to 400 mTorr, with 
an advanced feedback control of the current and a custom developed 12-bit ultra- fast imaging 
system with a timing resolution <1 msec. We have employed real-time 3D imaging acquired 
per laser pulse to monitor both specular and diffuse intensities. By recording the temporal and 
thermal evolution of the diffracted intensity during the growth combined with ex-situ atomic force 
microscopy (AFM), we were able to obtain exquisite details of the evolution of layer morphology. 

At first, we demonstrate the difference between the continuous low-frequency deposition (Figs. 
l(a, b)) and interrupted high-frequency (Figs. (Ic, d)) growth[7l |91 of the strongly correlated 
SLs. Figure 1(a) shows the result of a typical continuous deposition of ~4 uc thick layer of 



LNO with an imposed slow growth rate defined by the laser frequency of 3 Hz. As seen, the 
RHEED specular intensity (RSI) exhibits strong damping right after a few monolayers attesting to 
the rapidly increasing surface disorder. The diffuse intensity (inset in Fig. 1(a)) shows a series of 
pronounced peaks typical of electrons transmitted through microcrystalline 3D-like islands [fTOl . 
The large area 1/im x lyum AFM scans shown in Fig. 1(b) further corroborate the presence of 
pronounced 3D microstructures. Figures l(c, d) present the results obtained with the deposition 
rate increased up to 30 laser-pulses/sec and introduction of a prolonged time-delay (dwell time) 
up to 100 seconds between two consecutive unit cells. 

As immediately evident from Figs. l(c, d), such a LBL growth yields perfect coverage and an 
excellent morphology (see also auxiliary supplement[llj). We further performed detailed tempo- 
ral and temperature (T) dependent studies (from 450°C to 780°C) of the layer coverage during the 
growth of a sub-monolayer. The process of nucleation and condensation at high laser frequency re- 
vealed a non-trivial T-dependence. Specifically, the initial stage of deposition (area 1 in Fig. 2(a)) 
is T-independent, within the studied temperature interval, followed by two distinct non-linear T- 
dependences in the recovery stage 2 and the saturation stage 3. The origin of this interesting 
behavior can be connected to the specifics of the nucleation and condensation processes [fT2l [T3ll 
at the limit of high supersaturation and flux modulation. Under the conditions in use, the high 
nucleation rate results in a small critical nucleus radius [fT4l [TSl . ic ~1, which is on the order of 
a unit-cell. Secondly, we emphasize that a significant increase in the laser frequency pumps the 
ablated material into the vicinity of the substrate and thus further enlarges the supersaturation that 
scales with the mass density (for more data see auxiliary supplement[l IJ and Ref. l|T6l ). 

To corroborate the results deduced from the diffraction data, we have performed ex-situ AFM 
imaging of the surface morphology by interrupting the growth at ~0.3, 0.5 and 0.7 of layer cov- 
erage, 6. The obtained morphological progression is shown in Figure 2(b). As seen in Figure 
2(b-left), low layer coverage, ~0.36', (stage 1) is characterized by a large number of small-scale 
non-percolated islands which are characteristic of immobile clusters with a large density. At the 
half coverage step, 0.5^, (Fig. 2(b-center)) a characteristic wavy pattern of the deposited material 
has been observed, which is the result of islands coalescing to form a low-dimensional percolated 
structure. This structure can possibly be described in the framework of the diffusion-limited ag- 
gregation theory[14J. As the deposition continues (see Figure 2(b-right)), the growth approaches 
the layer completion, 0.76^, and looks morphologically imperfect implying that the amount of de- 
livered material is not sufficient to fill 'voids' in the under-layed layer. 



To gain further insight into the kinetics of growth, we have calculated the surface coverage, 
9, and roughness. A, in the framework of the bilayer model [fTOl [TTl [TSl . We begin with the 
general expression for the RSI I{t)Mi [13, /(t) = /o[l - 26*1 + 26*2]^ where 6*1 and 6*2 are the 
time-dependent coverages of the 1st and the 2nd atomic layers, respectively, and Iq is a scaling 
factor. In addition, the coverage is constrained to 9i + 62 = mt, where m is the deposition rate, 
which vanishes as the laser is switched off (i.e. m = at t > t^;). Next, we set mtE = 1 implying 
that 02, if present during the deposition, should vanish during the dwell time. The corresponding 



roughness can be calculated as A(t) = JOi{l — 61 — 62^ + 62(2 — 61 — 6*2)^. Under these con- 
ditions, the experimental RSI have been used to obtain 9i, 62 and A as a function of time and 
temperature. The results of the numerical treatment are presented in Figs. 3 (a, b). As anticipated 
for all the temperatures, the coverage of the 1st layer scales with the deposition time and is clearly 
T- independent during stage 1 (see Fig. 3(a)). On the other hand, 62 shows T-dependence only 
as the RSI approaches its lowest value at Im- In addition, the roughness, A shown in Fig. 3(b) 
increases rapidly from the start of deposition. Based on these observations, we infer that during 
stage 1 the kinetics-driven condensation is dominant over the substrate temperature and is crucial 
for the high-quality LBL growth. 

Starting from the point Im, the data obtained between 450°C and 780°C are markedly T- 
dependent. The SLs grown at around 700°C and above continue the same tendency found in 
stage 1, namely the continuing increase of 6*1 accompanied by the rapid decline of 62- In sharp 
contrast, for the samples grown at 450°C and SSO^C, 6*1 becomes progressively slower, and instead 
62 continues its strong progression. This behavior attests that the temperature-controlled diffusion 
becomes a key factor to maintain 2D growth and to suppress the development of 3D islands inher- 
ent to the colder deposition. The strong T-dependence of O2 also implies that for high temperatures 
the occasionally grown 2nd level is able to rapidly relax to the 1st level, whereas for colder deposi- 
tion this process is strongly suppressed resulting in the rapidly increasing roughness. Additionally, 
strong evidence for the new scaling law is deduced from an analysis of the T-dependence of the 
RSI after Im- Indeed, the T-independent scahng of stage 1 has now switched to a power law with 
exponent exponent (3 ^1/4 (further details in auxiliary supplement[jTIl). This marked difference 
can be further explicated from the observation that at Im, di ~0.5. This specific value is known to 
be critical in the 2D percolation problem[12J. At this point the new spatial arrangement appears 
and is characterized by a developing long-range order parameter as experimentally evidenced by 
the recovering RSI and characteristic morphology of percolated islands shown in Fig 2(b). 



Finally, upon completion of the ablation at t = tE (stage 3 in Fig. 3(a)), the temporal depen- 
dence of the coverage changes again. During this relaxation period, 62 starts decreasing, since the 
1st layer recovers at the expense of the 2nd layer. However, this recovery is strongly T-dependent. 
As seen, for the low temperature deposition 9i cannot recover in full, whereas at 700°C the re- 
covery is practically complete. During this phase the disorder configuration undergoes yet another 
change towards a new power law with (3 ^1. The experimental observation of the scaling specific 
to a particular mechanism of growth is important, because during each stage the system evolves by 
changing only a mean size of the surface disorder and without switching the disorder distribution 
function defined by the growth mode liT9ll20ll . 

In summary, we have developed the LBL growth of a novel lu LNO/luc LAO HS on STO 
substrate with excellent morphological quality. We have demonstrated the profound difference 
between the kinetics of continuous and interrupted growth of this HS and have investigated differ- 
ent stages of sub-monolayer nucleation and growth at the limit of high supersaturation and rapid 
flux modulation. These experimental findings are of fundamental and practical importance and 
should enable synthesis of unusual metastable material phases of novel ultra-thin complex oxide 
heterostructures. 
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FIG. 1: (Color onUne) (a) RSI oscillations of LNO/LAO SL (n?«6) grown at T = 730°C, Pq^ ~100 mTorr 
and with a laser power density ~2.2 J/cm^. Continuous deposition with a laser frequency of 3 Hz. (c) 
Interrupted growth with a frequency of 30 Hz and a dwell time ~100 sec. Insets: (a,c) RHEED patterns and 
(b,d) corresponding AFM images. 
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FIG. 2: (Color online) (a) T-dependence of the RSI during the unit cell growth. The red arrows indicate the 
start, ts and the end, Ie of the ablation. The lowest diffracted intensity is marked as tjv/- Red, green and 
blue circles correspond to ~0.3, 0.5 and 0.7 coverage, (b) AFM scans (1 x 1/xm^) obtained by interrupting 
the growth at the corresponding coverage. Since the evolution of disorder configuration during stage 1 is 
independent of temperature, the laminae-like pattern is roughly representative of the same surface seen by 
RHEED. Insets show the magnified false colored phase-contrast images. 
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FIG. 3: (Color online) (a) Temperature evolution of surface coverage, 6*1^2 and (b) corresponding layer 
roughness, A. The stages of growth marked as 1-3 correspond to the same time-intervals in Fig. 2(a).Note, 
the difference between the minimum and maximum of RSI is normalized to unity. 
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Auxiliary information for "Sub- monolayer nucleation and growth of complex oxides 
heterostructures at high supersaturation and rapid flux modulation" 
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FIG. 1: Auxiliary Supplement: (A) The slow intensity decay along the extended specular (OOL) x-ray diffraction rod 
demonstrates the highly ordered nature of the surface and interface of the superlattice film structure. The interference fringes 
verify the strong coherent layered epitaxy of the film. (B) presents an in-film-plane (HK) reciprocal space map of the (222)fiini 
reflection showing the fully in-plane lattice constraint to that of the SrTiOs (001) substrate. 
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FIG. 2: Auxiliary Supplement: (Color online) A-B. Scaling behavior of the RHEED specular intensity (RSI) as a function 
of temperature showing three distinct stages 1-3 characterized by distinct power laws found during the unit cell formation. 
Solid lines represent a fit to ln(I /lo) ex T^ . 
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FIG. 3: Auxiliary Supplement: (Color online) Evolution of the activation energy, Ta during the relaxation stage 3 with 
time. 



ON TIME DEPENDENCE OF THE RSI. 



Even though a theoretical description of the nanoscale growth at the supersaturation (SS) limit is beyond modern 
growth theories, the time evolution of the layer coverage can be evaluated analytically. Indeed, we can detail the kinet- 
ics of the 2nd layer coverage as follows, d92{t)/dt = \6i(t) — 62/T. Here the 1st term implies that the coverage of the 
2nd layer scales with the available surface of the 1st layer, because the 2nd layer cannot develop in the vacuum. The 
2nd term describes the relaxation of the 2nd layer. Correspondingly, this model involves two phenomenological param- 
eters, A and t. The model has an analytical solution, which for t < ts reads as 62 (t) = {\m/a)t—{\m/a'^) [1 — e^"*]. 



here a = A + 1/t, with the 1st non-zero term 92{t < t^) ex: t 
the exponentially decreasing coverage 62 (t) 



g(l) 



(t 



E e 



a(t-tE) 



During the dwell time (m = 0) the model yields 
As before, the corresponding 6i{t) shows that for 
t<tE 9i{t) increases linearly while 02(i) increases quadratically with time. After t = tE, di{t) continues growing while 
02 (t) rapidly declines. These model predictions are in accord with the observations. 



ON CRITICAL ROLE OF THE LASER FREQUENCY ON GROWTH KINETICS 

Here we discuss the critical role of the laser frequency on the kinetics of the early stage of growth beyond the 
nucleation timescale. As seen in Fig.A4 there are three distinct time scales tp, 1/ f and the time between deposition 
and aggregation of particles, ir - all three determining the kinetics of the growth between two consecutive pulses. 
First we note that a significant increase in the laser frequency pumps the ablated material into the vicinity of the 
substrate. This enlarges the SS which scales with the mass density. After that, a large number of monomers rapidly 
appear during the "flux-on" phase, tp, which for our specific experimental conditions lasts ^ 10~^-10~^sec[l].The 
surface coverage continues to evolve, and continue to evolve, until the inonomers reach a steady-state concentration 
due to adatom-adatom aggregation. 

Under the studied PLD conditions, this diffusion driven aggregation is characterized by the typical relaxation time 
tr (x 1/DN ~ 10^^-lO^^sec, where the island density, N and the diffusion constant, D are estimated from the fits of 
the RSI and data from Ref.s[2, 3]. For the low pulse-rate growth characterized by tp <^ tr, the monomers rapidly and 
within a single 'flux-on' period reach the steady-state concentration with no further aggregation occurring between 
two consecutive pulses. Upon increasing / beyond lOHz (i.e. 1// ~ 3 x 10^^-lO^^sec), the growth crosses into a 
new regime where tr ~ 1// (see Fig.A4). In this case, the transitional kinetics imposed by the rapid laser pulsing 
result in a monomer concentration which does not fully decay during the time between pulses. As a result, it takes 
several successive pulses to achieve the steady state. A careful scaling analysis of this transitional kinetics regime 
for the island size i=l also revealed that saturation of island density is independent of diffusion[4, 5]. As shown in 
Fig. 2(A), this is in excellent agreement with the observed T-independent behavior of the RSI during the stage 1. 
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FIG. 4: Auxiliary Supplement: (Color online) Timing diagram showing the inodulated flux (red bars) and relaxation of the 
monomer density (broken lines) (adapted from Ref. [4]). 

Notice, upon further increase of the pulse rate (i.e. tp <^ 1//) the growing system will 'see' only an average flux, _Fav 
and now fall into the regime resembling a continuous MBE-like growth[6]. Based on those observations, we conclude 
that pulse-rate plays the role of a new powerful kinetic 'handle' analogous to temperature albeit without activated 
behavior of the latter. 
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